Polyominoes with nearly convex columns: 
A semi-directed model 

Svjetlan Feretic * 
Faculty of Civil Engineering, University of Rijeka, 
Viktora Cara Emina 5, 51000 Rijeka, Croatia 

November 23, 2010 



Abstract 

Column-convex polyominoes are by now a well-explored model. So far, however, no 
attention has been given to polyominoes whose columns can have either one or two con- 
nected components. This little known kind of polyominoes seems not to be manageable 
as a whole. To obtain solvable models, one needs to introduce some restrictions. This 
paper is focused on polyominoes with hexagonal cells. The restrictions just mentioned 
are semi-directedness and an upper bound on the size of the gap within a column. The 
solvable models so obtained have rational area generating functions, as column-convex 
polyominoes do. However, the growth constants of the new models are 4.1149 . . . and 
more, whereas the growth constant of column-convex polyominoes is 3.8631 .... 
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1 Introduction 



The enumeration of polyominoes is a topic of great interest to chemists, physicists and 
combinatorialists alike. In chemical terms, any polyomino (with hexagonal cells) is a possi- 
ble benzenoid hydrocarbon. In physics, determining the number of n-celled polyominoes is 
related to the study of two-dimensional percolation phenomena. In combinatorics, polyomi- 
noes are of interest in their own right because several polyomino models have good-looking 
exact solutions. 

Let a n be the number of n-celled polyominoes. An exact formula for a n seems unlikely 
to ever be found. However, there exist notable results on the quantity c = lim„_ i . oc ^fa^. 
(This quantity is called the growth constant.) In the case of polyominoes with hexagonal 
cells, Voge and Guttmann [12] gave a rigorous proof that the growth constant c exists 
and satisfies the inequality 4.8049 < c < 5.9047 . In the same paper, it is estimated that 
c = 5.1831478(17). 

In this paper, we are not going to improve the above-stated bounds on c. Instead, we are 
going to revisit column- convex polyominoes. The first to study this now-familiar model was 
Temperley [11] in 1956. Whether with square cells or with hexagonal cells, column-convex 
polyominoes have a rational area generating function. When cells are hexagons, the growth 
constant of column-convex polyominoes is 3.8631 . . . (Klarner [10 ). This growth constant 
remained a record (among polyomino models having reasonably simple exact solutions) 
until 1982, when Dhar, Phani and Barma [3] invented directed site animals. Directed site 
animals with step-set A = {(1, 0), (0, 1)} can be viewed as a family of polyominoes with 
square cells, while directed site animals with step-set B = {(1, 0), (0, 1), (1, 1)} can be 
viewed as a family of polyominoes with hexagonal cells. If the step-set is B, the growth 
constant of directed site animals is exactly 4. Incidentally, with either of the step-sets A 
and B, the area generating function of directed site animals is algebraic (which is rather 
surprising) and satisfies a simple quadratic equation. Later on, in 2002, Bousquet-Melou 
and Rechnitzer [2 introduced stacked directed animals and multi- directed animals. Those 
two models substantially generalize directed site animals. Whether the cells are squares or 
hexagons, the area generating function of stacked directed animals is degree-two algebraic, 
and the area generating function of multi-directed animals is not D-finite. When cells are 
hexagons, the growth constant of stacked directed animals is exactly 4.5, and the growth 
constant of multi-directed animals is 4.5878 .... 

Although descended from directed animals, multi-directed animals are also a superset of 
column-convex polyominoes. (To be precise, this holds when the cells are hexagons. It is 
not quite clear whether multi-directed animals with square cells are a superset of column- 
convex polyominoes with square cells.) Besides multi-directed animals, there exist two 
other generalizations of column-convex polyominoes. Those two generalizations are called 
m-convex polygons [S] and prudent polygons [5]. So far, however, m-convex polygons and 
prudent polygons have not been enumerated by area; they have only been enumerated by 
perimeter. 

The aim of this paper is to define a model which (a) generalizes column-convex polyomi- 
noes, (b) possesses a reasonably simple area generating function, and (c) possesses a high 
growth constant. In view of the facts stated above, we shall have to compete with just one 
pre-existing model, namely with multi-directed animals. 

In this paper, we actually introduce a sequence of generalizations of hexagonal- celled 
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column-convex polyominoes. Namely, we define level to cheesy polyominoes (to = 1,2,3,...). 
At every level, our new model has a rational area generating function. However, as level 
increases, those rational generating functions quickly gain in size. In computing generating 
functions, we go up to level 3. The computations are done by using Bousquet-Melou's [T] 
and Svrtan's [7] "turbo" version of the Temperley method [TT]. The growth constants of 
level 1, level 2 and level 3 cheesy polyominoes turn out to be 4.1149 . . . , 4.2318 . . . and 
4.2886 . . ., respectively. 

Thus, the growth constants of cheesy polyominoes are not as high as 4.5878 . . ., the 
growth constant of multi-directed animals. However, we expect to obtain greater growth 
constants (4.5 or more) in our next two papers. Namely, we have in view two generaliza- 
tions of level to cheesy polyominoes. The names of those two generalizations are level to 
polyominoes with cheesy blocks and level to column- subconv ex polyominoes. Incidentally, for 
every to € N, the area generating function of level m polyominoes with cheesy blocks is 
rational, whereas the area generating function of level to column-subconvex polyominoes is 
probably not D-finite. We are also planning to study similar generalizations of square-celled 
column-convex polyominoes. Let us mention, however, that our new hexagonal-celled mod- 
els behave somewhat better than their square-celled counterparts. Solving a square-celled 
level 1 model requires almost as much effort as solving a hexagonal-celled level 2 model. 

2 Definitions and conventions 

There are three regular tilings of the Euclidean plane, namely the triangular tiling, the 
square tiling, and the hexagonal tiling. We adopt the convention that every square or 
hexagonal tile has two horizontal edges. In a regular tiling, a tile is often referred to as a 
cell. A plane figure P is a polyomino if P is a union of finitely many cells and the interior 
of P is connected. See Figure 1. Observe that, if a union of hexagonal cells is connected, 
then it possesses a connected interior as well. 

Let P and Q be two polyominoes. We consider P and Q to be equal if and only if there 
exists a translation / such that f(P) = Q. 

From now on, we concentrate on the hexagonal tiling. When we write "a polyomino" , 
we actually mean "a hexagonal-celled polyomino" . 

Given a polyomino P, it is useful to partition the cells of P according to their horizontal 
projection. Each block of that partition is a column of P. Note that a column of a polyomino 
is not necessarily a connected set. An example of this is the highlighted column in Figure 1. 
On the other hand, it may happen that every column of a polyomino P is a connected set. 
In this case, the polyomino P is a column-convex polyomino. See Figure 2. 

Let a and b be two adjacent columns of a polyomino P. Let the column b be the right 
neighbour of the column a. It is certain that column b has at least one edge in common 
with column a; otherwise P could not be a connected set. However, there is no guarantee 
that every connected component of b has at least one edge in common with a. For example, 
in Figure 1, the upper component of the highlighted column has no edge in common with 
the previous column. 

Suppose that P is such a polyomino that the first (i.e., leftmost) column of P has no 
gap and that, in every pair of adjacent columns of P, every connected component of the 
right column has at least one edge in common with the left column. Then we say that P is 
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Figure 1: A hexagonal-celled polyomino. 




Figure 3: A level one cheesy polyomino. 



a rightward- semi- directed polyomino. 

A polyomino P is a level m cheesy polyomino if the following holds: 

• P is a rightward-semi-directed polyomino, 

• every column of P has at most two connected components, 

• if a column of P has two connected components, then the gap between the components 
consists of at most m cells. 

See Figure 3. 

Level one cheesy polyominoes are a subset of level two cheesy polyominoes, level two 
cheesy polyominoes are a subset of level three cheesy polyominoes, and so on. As m tends 
to infinity, the set of level m cheesy polyominoes tends (in a certain sense) to the set of all 
polyominoes which are rightward-semi-directed and are made up of columns with at most 
two connected components. 

The name "cheesy polyominoes" is intended to suggest that these polyominoes can have 
internal holes. At this point, it might be objected that there exists another model, called 
directed animals, in which holes occur more freely than in cheesy polyominoes. A column of 
a directed animal can have any number of holes, and the sizes of those holes are not subject 
to any limitations. However, the name of our model has 11 years of traditiorQ. So, our 
model will retain the name "cheesy polyominoes" despite the fact that directed animals are 
arguably "cheesier". 

If a polyomino P is made up of n cells, we say that the area of P is n. 

lr The contents of this paper exist since 1999. In that year, I defined cheesy polyominoes, explored them, 
and presented them at Mathematical Colloquium in Osijek [4], as well as at MATH/CHEM/COMP Course 
& Conference in Dubrovnik [5]. 
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Let a be a column of a polyomino P. By the height of a we mean the number of those 
cells which make up a plus the number of those (zero or more) cells which make up the gaps 
of a. For example, in Figure 1, the highlighted column has height 7, and the next column 
to the left has height 4. 

Let R be a set of polyominoes. By the area generating function of R we mean the formal 
sum 

^ ^ qCirea of P 
PER 

By the area and last column generating function of R we mean the formal sum 

^ ^ qCirea of P jthe height of the last column of P 
PER 



3 Cheesy polyominoes vs. multi-directed animals 

This section has been removed because arXiv warned us that the submission exceeds size 
limits. The removed section contained definitions and some pictures of the directed classes 
mentioned in the introduction. From those definitions and pictures, it was clear that cheesy 
polyominoes are very different from the said directed models. 



4 Column-convex polyominoes 

The area generating function for column-convex polyominoes is known since 1967 (Klarner 
[10]). However, rederiving that formula here will add to the completeness of this paper. 

Let A = A(q,t) be the area and last column generating function for column-convex 
polyominoes. Let A\ = A(q, 1) and B\ = ^(q, 1). 

Let S be the set of all column-convex polyominoes. 

We write S a for the set of column-convex polyominoes which have only one column. For 
P G S \ S a , we define the pivot cell of P to be the lower right neighbour of the lowest cell 
of the second last column of P. See Figure 9. Observe that the pivot cell of a polyomino 
P (E S \ S a is not necessarily contained in P. 

Let 



Sp — {P £ S\ S a : P contains its pivot cell} and 
5 7 = {P € S \ S a : P does not contain its pivot cell}. 

The sets S a , Sp and i5 7 form a partition of 5*. We write A a , Ap and A-y to denote the 
parts of the scries A that come from the sets S a , Sp and S 7 , respectively. 
It is obvious that 

A a =qt + (qt) 2 + (qtf + ... = -^-. (1) 



G 



The gray cell 
is the pivot cell. 



Figure 9: The pivot cell. 



If a polyomino P lies in Sp, then the last column of P is made up of the pivot cell, of 
i e {0, 1, 2, 3, . . . } cells lying below the pivot cell, and of j e {0, 1, 2, 3, . . . } cells lying 
above the pivot cell. Hence, 



Ax ■ qt 



1=0 



3=0 



qt 



(i - q ty 



(2) 



See Figure 10. 

For n € {1, 2, 3, . . .}, let denote the set of column-convex polyominoes whose 

last column consists of n cells. Let A^ be the part of A\ that comes from the set S (n \ 
Now, every element of 5 7 can be produced in three steps. Step one: We choose a number 
n e {1, 2, 3, . . .} and a polyomino P e S^ n \ Step two: In the last column of P, we choose 
a cell c. Step three: After choosing a number i e {1, 2, 3, . . .}, we place a new column of 
height i so that the lowest cell of the new column is the upper right neighbour of the cell c. 
See Figure 11. Thus, 



A, = E^ ) --E(9») i = E n -^ 1 -i 



qt 



qt 



qt 



1 - qt 

qt 

1 - gt 
qt 



[A? + 2A& + 3A^ + ... 



4 1) + 24 2) t+34 3) t 2 + .. 



with t—1 



dA(g,t) 
dt 



qt 



with t—1 



1 - qt 



1 - qt 

Since A = A a + Ap + A 1 , Eqs. (1), (2) and (3) imply that 

qt qt qt 



■Bx. 



A 



+ 



l-qt (1 - qt)* 



■Ax 



1 - qt 



Bx. 



(3) 



(4) 
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Figure 10: The last two columns of an element of Sp. 
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Setting t = 1, from Eq. (4) we obtain 



M = 



■A 1 



■Si 



1-9 (1-9) 2 ~~* 1-9 
Differentiating Eq. (4) with respect to t and then setting t = 1, we obtain 

9 , 2g 2 



Si 



1-9 (1-9) 2 



(1-9)2 (l-g)S 



1-9 (1-9) 2 



We now have a system of two linear equations in two unknowns, A\ and B\. By solving 
the system, we get the following proposition. 

Proposition 1 The area generating function for column- convex polyominoes is given by 

9'(l-9) 3 



A 1 



1 - 6q + 10q 2 - 7q 3 



Being a quartic polynomial, the denominator of A\ has four roots. Those roots are 
n = 0.2588 r 2 = 0.7066 ... - 0.4750 . . . • i, r 3 = 0.7066 . . . + 0.4750 . . . • i and r 4 = 
5.3278 .... The root with smallest absolute value is n = 0.2588 . . ., and A- is equal to 
3.8631 .... Therefore, the coefficient of q n in A\ (denoted [g™]Ai) has the asymptotic 
behaviour [g n ]Ai <~ c x 3.8631 . . .", where c is a constant. To find the value of c, we 
decompose A\ into partial fractions. The partial fraction involving n turns out to be 
In the Taylor series expansion of — °n t^ss > l nc coefficient of q n is equal to 



0.0487.. 



qr-0.2588... 

0.1884... x 3.8631 



g-0.2588. 

Thus, c = 0.1884 .... We have got the following result. 



Corollary 1 The number of n- celled column-convex polyominoes \q n \A\ has the asymptotic 
behaviour 

[q n ]A 1 - 0.188419 ... x 3.863130 



Thus, the growth constant of column-convex polyominoes is 3.8631 .... 



5 Level one cheesy polyominoes 

Counting level one cheesy polyominoes by area is nearly as easy as counting column-convex 
polyominoes by area. 

Let C = C(q,t) be the area and last column generating function for level one cheesy 
polyominoes. Let C\ — C(q, 1) and D\ — ^r{q, 1)- 

Let T be the set of all level one cheesy polyominoes. We write T a for the set of level one 
cheesy polyominoes which have only one column. For P e T\T a , we define the pivot cell 
of P to be the lower right neighbour of the lowest cell of the second last column of P. As 
with column-convex polyominoes, the pivot cell of a polyomino P £ T\T a is not necessarily 
contained in P. Let 
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Tp = {P £T\T a : the last column of P has no hole, 

and the pivot cell of P is contained in P}, 
Tj = {P £ T\T a : the last column of P has no hole, 

and the pivot cell of P is not contained in P}, and 
T s = {P £ T \ T a : the last column of P has a hole}. 

The sets T a , Tp, T 1 and Ts form a partition of T. We write C a , Cp, C 7 and Cs to denote 
the parts of the series C that come from the sets T a , Tp, T 7 and Tg, respectively. 
Similarly as with column-convex polyominoes, we have 



C' 



qt + (qt) 2 + (qt) 3 + ... 



qt 



1 - qt 



and 



Cp = 



qt 



■Ci- 



(5) 



(6) 



(l-qtf 

Consider the following situation. A cheesy polyomino P ends with a column I. We are 
creating a new column to the right of /, and the result should be an element of T~. Then, 
whether or not the column / has a hole, we can put the lowest cell of the new column in 
exactly m places, where m is the height of I. See Figure 12. Hence 



qt 



■ D 1 . 



(7) 



1 - qt 

Let us move on to another situation. A cheesy polyomino P ends with a column J. We 
are creating a new column to the right of J, and the result should be an element of Ts. 
Then, whether or not the column J has a hole, we can put the hole of the new column in 
exactly n—\ places, where n is the height of J. Sec Figure 13. 

For r £ {1, 2,3,.. .}, let T< r ) denote the set of level one cheesy polyominoes whose last 
column has height r. Let C[ r ^ be the part of C\ that comes from the set T^ r \ Every element 
of Ts can be produced in four steps. Step one: We choose a number r £ {1, 2, 3, . . .} and 
a polyomino P £ T^ r \ Step two: In the last column of P, we choose two adjacent cells. 
(This can be done in r — 1 ways.) Step three: We choose two numbers, i and j, from the 
set {1, 2, 3, . . .}. Step four: We create a two-component new column in which the upper 
component has i cells, the lower component has j cells, and the hole is the common right 
neighbour of the two cells chosen in Step two. Thus, 



r=l 



,i=l 



r— 1 ^ 



2 + 3 



(r) qt 



(l-qtf 



i=i 
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jure 12: The last two columns of two elements of T, 



7- 



12 




Figure 13: The last two columns of two elements of T$. 
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2.3 



qt 



(1 - qtf 
q 2 t 3 



cf ) +2Cf ) +3C 1 (4) + ... 
{i _ qt)2 {[c^ + 2Cfk + 3Cf¥ + ... 

c[ 1) t + c[ 2) t 2 + c { 1 3 h :i + ..] t } 



unih i— 1 



2+3 



(1 - gi) 2 
g 2 t 3 



dC{q,t) 



dt 

■(Di-Ci). 



with t— 

[ c (q,t)] with *=i 



(1 - gi) 2 

Since C = C Q + + C 7 + C 5 , Eqs. (5), (6), (7) and (8) imply that 



C = 



qt 



+ 



qt 



Ci + 



1-gi (1-gt) 2 ~* ' 1-gf 
Setting i = 1, from Eq. (9) we obtain 



+ 



Ci + 



£>! + 



£>! + 



g 2 t 3 



(£>i - Ci). 



(1 - gi) 2 

• - Ci). 



(8) 



(9) 



1-9 ' (1-9) 2 ^ ' 1-9 ~* ' (I-?) 2 
Differentiating Eq. (9) with respect to t and then setting t = 1, we obtain 



+ 



1-9 (l-9) : 



+ 



+ 



9 



+ ■ 



2g 2 



(1-g) 2 (1-g) 3 
3g 2 



1-9 (1-9) 2 



■Ci 
2q 3 



(1-g) 2 



(Di - Ci). 



Thus, things are similar as with column-convex polyominoes. We have a system of 
two linear equations in two unknowns, C\ and D\. By solving the system, we get this 
proposition. 

Proposition 2 The area generating function for level one cheesy polyominoes is given by 

_ g(l-3g + g 2 ) 
1 1 - 6g + 8q 2 - q 3 ' 

The roots of the denominator of A\ are n = 0.2430 . . ., r 2 = 0.5727 . . . and r 3 = 
7.1842 .... (All the three roots are real.) The root with smallest absolute value is n = 
0.2430..., and ^- is equal to 4.1149... . By decomposing C\ into partial fractions and 
expanding the partial fractions into Taylor series, we find out the following fact. 

Corollary 2 The number ofn-celled level one cheesy polyominoes [q n ]C\ has the asymptotic 
behaviour 

[<7 ri ]Ci - 0.144176 ... x 4.114907 

Thus, the growth constant of level one cheesy polyominoes is 4.1149 .... We observe 
a considerable increase with respect to 3.8631 . . ., the growth constant of column-convex 
polyominoes. 
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6 Level two cheesy polyominoes 



In this enumeration, if the last column of a polyomino has two connected components, we 
sometimes need to record not only the overall height of the last column, but also the height 
of the last column's upper component and the height of the last column's lower component. 
Hence, in addition to the "old" variables q and t, we introduce two new variables, u and v. 
As before, the exponent of q is the area and the exponent of t is the overall height of the 
last column^). The exponent of u is the height of the upper component of the last column, 
and the exponent of v is the height of the lower component of the last column. 

The two main generating functions in this enumeration are E = E(q,t) and G = 
G(q,t,u,v). Those generating functions are used for the following purposes: 

• E is a generating function for level two cheesy polyominoes whose last column either 
has no hole or has a one-celled hole, 

• G is a generating function for level two cheesy polyominoes whose last column has a 
two-celled hole. 

Let E x = E(q,l), F = $f(g,0), F l = ^(q,l), G 1 = G(q, 1,1,1), H x = ™(q, 1,1,1), 
/„ = 1,0, 1), and J = 1, 1, 0). 

Let U be the set of those level two cheesy polyominoes whose last column either has no 
hole or has a one-celled hole. Let V be the set of those level two cheesy polyominoes whose 
last column has a two-celled hole. For P G U U V, we dehne the body of P to be all of P, 
except the rightmost column of P. 

We are now going to partition the set U into seven subsets and the set V into two subsets. 

Let 

U a = {P € U : P has only one column}. 

Up = {P E U\U a : the body of P lies in U, the last column of P has no hole, 

and the pivot cell of P is contained in P}, 
U~f = {P € U \U a ■ the body of P lies in U, the last column of P has no hole, 

and the pivot cell of P is not contained in P}, 
Us = {P G U \U a : the body of P lies in U, and the last column of P has a hole}, 
U e = {P € U \ U a ■ the body of P lies in V, the last column of P has no hole, 

and the pivot cell of P is contained in P}, 
U{ = {P € U \U a : the body of P lies in V, the last column of P has no hole, 

and the pivot cell of P is not contained in P}, and 
U, q = {P G U \U a : the body of P lies in V, and the last column of P has a hole}. 

Let 

2 Recall what do we mean by the height of a column: in Figure 1, the highlighted column has height 7, 
and the next column to the left has height 4. 
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V a = {P e V : the body of P lies in U}, and 
Vp = {P e V : the body of P lies in V}. 

It is clear that the sets U a , Up,..., U v form a partition of U, and that the sets V a and 
Vp form a partition of V. We shall write E a , Ep, . . . , E„ for the parts of the series -E that 
come from the sets U a , Up, ... , ET^, respectively. Also, we shall write G a and Gp for the 
parts of the series G that come from the sets V a and Vp , respectively. 

As in Section 5, we have 

E a = (10) 
Es = J^-iF.-E^ (13) 

E * ^ (rV Gl - (14) 

The functional equation for is more interesting. Let P be an clement of V, let I be 
the last column of P, and let m be the height of I. Suppose that we are creating a new 
column to the right of I, and that the result should be an element of U^. To the right of I 
there is one cell (say c) which shares an edge with each of the two cells forming the hole of 
/. If we choose c as the lowest cell of the new column, then the new column will have to 
be at least two cells high. Otherwise P U (the new column) will not be a polyomino. Sec 
Figure 14. In addition to c, there are m — 1 other choices for the lowest cell of the new 
column. For each of these m — 1 choices, P U (the new column) is a polyomino regardless 
of how many cells the new column has. Altogether, we have 

We proceed to E v . Let P be an element of V, let J be the last column of P, and let 
n be the height of J. Suppose that we are creating a new column to the right of J, and 
that the result should be an element of U v . One of the choices for the hole of the new 
column is the upper right neighbour of the top cell of the lower component of J. For this 
choice, the upper component of the new column has to have at least two cells. Otherwise 
P U (the new column) is not a polyomino. See Figure 15. We can also choose the hole of 
the new column as the lower right neighbour of the bottom cell of the upper component 
of J. Then, in order for P U (the new column) to be a polyomino, the lower component of 
the new column has to have at least two cells. In addition to the two ways just considered, 
there exist n — 3 other ways to choose the hole of the new column. For each of those n — 3 
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(a) (b) 

Figure 14: (a) For this choice of the lowest cell, the new column must be at least two cells 
high, (b) For this choice of the lowest cell, the new column does not have to be at least two 
cells high. 
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Figure 15: (a) If the one-celled hole is in this position, the upper component of the new 
column must have at least two cells, (b) If the one-celled hole is in this position, the upper 
and lower components of the new column can be of any sizes. 



ways, P U (the new column) is a polyomino regardless of the sizes of the upper and lower 
components of the new column. Thus, we have 

2o 3 i 4 a 2 t 3 

E *=(dw- Gl+ (^w- {Hl - 3Gl) - (16) 

Of course, the series E a , Ep,. . ., E v sum up to E. Therefore, the summation of Eqs. 
(10)-(16) gives 

F qt l qt F l qt F i (F F) 
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Figure 16: The last two columns of two elements of V a . 



+ 



qt 



■Gi + 




• Gi + 



qt 



■ (#1 - GO 



+ 



(l-qtf 
2q 3 t 4 



■G t + 



l-qt' 
qH 3 



■(^ -3d). 



1 - qt 



(17) 



(1 - q ty 



(1 - qty 



We also need to establish a functional equation for G. Let P be an element of U and let 
c be a column with a two-celled hole. Suppose that we want to glue c to P so that PL) c lies 
in V a , and so that P and c are the body and the last column of P U c, respectively. In how 
many ways P and c can be glued together? In principle, the number of ways is (the height 
of the last column of P) minus two. See Figure 16. However, if P ends with a one-celled 
column, then we can glue c to P in zero ways, and not in minus one ways. Thus, we have 



In the case of Gp, it is convenient to use overcounting. That is, we are going to "mis- 
takenly" count too much, and then subtract the parts which do not belong. Let P be an 



G a — 




■ {F 1 -2E 1 +F ). 



(18) 



(1 - qtu)(l - qtv) 
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element of V and let c be a column with a two-celled hole. Suppose that we want to glue 
c to P so that P U c lies in Vp , and so that P and c are the body and the last column of 
P U c, respectively. In how many ways P and c can be glued together? First, there are (the 
height of the last column of P) minus two ways to satisfy these two necessary conditions: 

• the bottom cell of the upper component of c is either identical with or lies lower than 
the upper right neighbour of the top cell of the last column of P, and 



the top cell of the lower component of c is either identical with or lies higher than the 
lower right neighbour of the bottom cell of the last column of P. 



See Figure 17. 

So, if there were no special cases, then Gp would be equal to 



q t uv 



(tfi-2Gi). (19) 



(1 - qtu){\ - qtv) 
However, special cases do exist. There are two of them: 

1. The upper component of the last column of P e V has at least two cells and the lower 
component of the two-component column c has just one cell. 

2. The lower component of the last column of P € V has at least two cells and the upper 
component of the two-component column c has just one cell. 

In case 1, it is (so to speak) dangerous to glue c to P in such a way that the one-celled 
lower component of c becomes a common neighbour of the two cells which form the hole of 
the last column of P. This dangerous operation produces an object which is not a polyomino 
and hence does not lie in Vp. 

In case 2, it is dangerous to glue c to P in such a way that the one-celled upper component 
of c becomes a common neighbour of the two cells which form the hole of the last column of 
P. Again, the dangerous operation produces an object which is not a polyomino and hence 
does not lie in Vp. 

Now, Eq. (19) is actually a generating function for the union of Vp with the set of 
objects produced by the two dangerous operations. The generating function for the objects 
produced by the first dangerous operation is 

The generating function for the objects produced by the second dangerous operation is 

PD 

Subtracting Eqs. (20) and (21) from Eq. (19), we obtain 

qH 4 uv q 2 t 4 uv q 2 t 4 uv 

Gf> ~ (1 - qtu)(l - qtv) ■ (Hl - 2Gl) - l^Tu ■ (Gl l0) l^qTv ' (Gl " ^ (22) 
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Figure 17: The last two columns of five elements of Vp. 
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Since G — G a + Gp, Eqs. (18) and (22) imply that 



G = 




■ (Fi — 2E\ + F + Hx — 2Gi) 



-^■(d-Io)- 
1 - qtu 



1 — qtv 



q 2 t 4 uv 



■ (Gi - Jo)- 



(23) 



Using the computer algebra system Maple, from Eqs. (17) and (23) we obtained a system 
of seven linear equations in seven unknowns: E\, F , F\, G\, H\, I and Jo- Rather than 
write down these seven equations (some of which are a bit cumbersome) , here below we give 
a list of recipes. Recipe no. k (k = 1, 2, . . . , 7) tells how to obtain the kth equation of the 
linear system. 

1. In Eq. (17), set t = 1. 

2. Differentiate Eq. (17) with respect to t and then set t = 0. 

3. Differentiate Eq. (17) with respect to t and then set t = 1. 

4. In Eq. (23), set t = u = v = 1. 

5. Differentiate Eq. (23) with respect to t and then set t = u = v = 1. 

6. Differentiate Eq. (23) with respect to u. Then set t = v = 1 and w = 0. 

7. Differentiate Eq. (23) with respect to v. Then set £ = u = 1 and w = 0. 

The computer algebra quickly solved the linear system and then summed the generating 
functions E\ and G\. The result can be seen in the following proposition. 

Proposition 3 The area generating function for level two cheesy polyominoes is given by 



Corollary 3 The number of n-celled level two cheesy polyominoes has the asymptotic be- 



Thus, the growth constant of level two cheesy polyominoes is 4.2318 .... 

7 Level three cheesy polyominoes 

The enumeration of level three cheesy polyominoes is similar to the enumeration of level two 
cheesy polyominoes. However, there are still more cases than before. In the enumeration of 
level two cheesy polyominoes, the total number of cases was 9 (because we partitioned the 
set U into 7 subsets and the set V into 2 subsets). At level three, the total number of cases 
is 16. We deem it reasonable to skip those 16 cases and only state the final result. 



K = 



q-{l-&q + llq 2 - 6q 3 - q 4 - 3q 6 + 5q 7 + 4q s - 3q 9 - 3q w ) 



l-9q+ 27q 2 - 3lq 3 + 8q A + 4<f - 2q e + 16q 7 - 5q 8 - 16q 9 - 2q w + 5q 



li 



haviour 



[q n ]K - 0.121042 ... x 4.231836 
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Proposition 4 

where 



The area generating function for level three cheesy polyominoes is given by 

L- M 
N , 



M = 9 • (1 - llg 



A9q 2 - 114g 3 



94<f + 5<f 

,13 



71q 7 - 143<? s 

,14 gi„15 



lA&q 1 

+ 176<? 9 - lhAq 10 + 100<? n + 24q 12 - 121g 13 + 90q 14 - 6lq 15 + 19q 16 
+ 58q 17 - 32q ls - 3lq 19 + 37q 20 + Uq 21 - A3q 22 - Aq 23 + 21g 24 
- q 25 - 5q 26 ) 



and 



N = l-Uq+ 80q 2 - 2A3q 3 + 423 g 4 - 413g 5 + 174g 6 + 106g 7 - 350g 8 
+ 533g 9 - 546g 10 + 427g n - 148g 12 - 261g 13 + 383g 14 - 253g 15 
+ 158g 16 + 57g 17 - 181g 18 + 10g 19 + llbq 20 - A9q 21 - 9Qq 22 
+ 93q 23 + A9q 24 - 5Aq 25 - 12q 2e + I2q 27 + q 28 . 

Corollary 4 The number of n- celled level three cheesy polyominoes has the asymptotic be- 
haviour 

[q n ]L ~ 0.108269 ... x 4.288630 



8 Taylor expansions and the limit value of the growth 
constants 

To see how many polyominoes of a given type have 1, 2, 3, . . . cells, we expanded the area 
generating functions into Taylor series. The results are shown in Table 1. 

The row "area = 4" of Table 1 reads 42, 43, 43, 43, 44. Indeed, Figure 18 shows the 
only two four-celled polyominoes which are not column-convex polyominoes. The polyomino 
on the left is a level m cheesy polyomino for every m £ N, and the polyomino on the right 
is not a level m cheesy polyomino for any m € N. 

Next: How do the growth constants behave when level tends to infinity? Our database 
is too small for making precise estimates. Anyway, we know that the growth constant 
of column-convex polyominoes is 3.863, while the growth constants of level one, level two 
and level three cheesy polyominoes are 4.115, 4.232 and 4.289, respectively. Computing 
the first differences, we get the numbers 4.115 - 3.863 = 0.252, 4.232 - 4.115 = 0.117, and 
4.289-4.232 = 0.057. Now, the sequence 0.252, 0.117, 0.057 is a little similar to a geometric 
sequence with common ratio \. Hence, the limit value of the growth constants of cheesy 
polyominoes might be about 4.232 + 2 • 0.057 = 4.346. 
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Column- 


Level 1 


Level 2 


Level 3 






convex 


cheesy 


cheesy 


cheesy 


All 




polyo- 


polyo- 


polyo- 


polyo- 


polyo- 


Area 


minoes 


minoes 


minoes 


minoes 


minoes 


1 


1 


1 


1 


1 


1 


2 


3 


3 


3 


3 


3 


3 


11 


11 


11 


11 


11 


4 


42 


43 


43 


43 


44 


5 


162 


173 


174 


174 


186 


6 


626 


705 


718 


719 


814 


7 


2419 


2889 


2996 


3012 


3652 


8 


9346 


11867 


12579 


12727 


16689 


9 


36106 


48795 


52996 


54067 


77359 


10 


139483 


200723 


223705 


230464 


362671 


11 


538841 


825845 


945324 


984477 


1716033 


12 


2081612 


3398081 


3997185 


4211222 


8182213 



Table 1: Here is how many polyominoes of a given type have 1, 2, . . . , 12 cells. 



Figure 18: The 43rd and 44th four-celled polyominoes. 
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9 Conclusion 



This paper is concerned with polyominoes on the hexagonal lattice. For every m € N, we 
have defined a set of polyominoes called level m cheesy polyominoes. A polyomino P is a 
level m cheesy polyomino if the following holds: 

1. P is a rightward-semi-directed polyomino, 

2. every column of P has at most two connected components, 

3. if a column of P has two connected components, then the gap between the components 
consists of at most m cells. 

Column-convex polyominoes are a subset of level one cheesy polyominoes and, for every 
m € N, level m cheesy polyominoes are a subset of level m + 1 cheesy polyominoes. 

For every m € N, level m cheesy polyominoes have a rational area generating func- 
tion. We have computed the area generating functions for levels one, two and three. At 
those three levels, the number of n-celled cheesy polyominoes is asymptotically equal to 
0.1441 ... x 4.1149 . . .", to 0.1210 . . . x 4.2318 . . .", and to 0.1082 . . . x 4.2886 . . .", respec- 
tively. For comparison, the number of n-celled column-convex polyominoes is asymptotically 
equal to 0.1884 . . . x 3.8631 . . .". The number of n-celled multi-directed animals behaves 
asymptotically as constant x 4.5878..." j2j. (At present, multi-directed animals are the 
largest exactly solved class of polyominoes. However, multi-directed animals are not a su- 
perset of level one cheesy polyominoes.) The number of all n-celled polyominoes behaves 
asymptotically as °- 2 ™- x 5. 1831 ..J 1 [12]. 

This work could be generalized in several ways. The requirement that polyominoes 
are rightward-semi-directed can be relaxed and even removed. However, when the above 
definition is reduced to requirements no. 2 and 3, the area generating function is not a 
rational function, but a complicated g-series. It is also possible not to require rightward-semi- 
directedness and, at the same time, allow two-component columns to have gaps of all sizes. 
The only remaining requirement is then the second one, "every column of P has at most 
two connected components" . The said requirement by itself defines an unsolvable model, 
but that model can be made solvable by introducing a new requirement. For example, if we 
forbid runs of two or more consecutive two-component columns, then the area generating 
function is again a complicated, but computable, g-series. 

As far as I can see, if requirement no. 2 is replaced by "every column of P has at most 
three connected components", the resulting model is still solvable, but if requirement no. 2 
is replaced by "every column of P has at most four connected components" , the resulting 
model is unsolvable. 

I think that, already at level one, cheesy polyominoes cannot be enumerated by perime- 
ter. Namely, the perimeter generating function has zero radius of convergence, as can be 
proved by adapting an argument given by Tony Guttmann in Section 9 of [6]. 
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